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Abstract 

Capillary jetting of a fluid dispersed into another immiscible phase is usually limited by a critical 
Capillary number, a function of the Reynolds number and the fluid properties ratios. Critical 
conditions are set when the minimum spreading velocity of small perturbations u* along the jet 
(marginal stability velocity) is zero. Here we identify and describe parametrical regions of high 
technological relevance, where u* > and the jet flow is always supercritical independently of the 
dispersed liquid flow rate: within these relatively broad regions, the jet does not undergo the usual 
dripping-jetting transition, so that either the jet can be made arbitrarily thin (yielding droplets of 
any imaginably small size), or the issued flow rate can be made arbitrarily small. In this work, we 
provide illustrative analytical studies of asymptotic cases for both negligible and dominant inertia 
forces. In this latter case, requiring a non-zero jet surface velocity, axisymmetric perturbation 
waves "surf" downstream for all given wave numbers while the liquid bulk can remain static. In 
the former case (implying small Reynolds flow) we found that the jet profile small slope is limited 
by a critical value; different published experiments support our predictions. 
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I. INTRODUCTORY REMARKS 



The quest for the conditions under which a given stream of fluid 1 can be dispersed as 
very small, homogeneously sized droplets into another immiscible fluid 2 is an old endeavor. 
Steady capillary jetting produces droplets of any desired diameter at a controllable rate 
through Rayleigh-Plateau instability, and thus is the preferred choice in many applications. 
Jetting cau be suited by a divers.ty of energy sources, from p.aiu pressure Q or elec- 
2j (or their combination |3|, |4|), to chemical potential [5j or even thermal 



ace under concentrated pho- 
6| or even under ultrasound 



trostatic suction 

gradients. Very recently, jetting has also been shown to take p 
ton irradiation (laser) when surface tension is extremely low 
irradiation by the same group. Capillary jetting from a fluid source gives rise to droplets 
smaller than dripping (a phenomenon where drops are individually issued from the source at 
a certain frequency) under the same operating conditions. Consequently, a significant effort 
las been lately devoted to map the transition from jetting to dripping from a fluid source 
3, S [2I, [lQ|, in the search for extended jetting conditions down to the smallest possible jet di- 
ameter. Two key dimensionless numbers gauge the role of inertia and viscous forces relative 
to surface tension, namely Weber and Capillary numbers We = piU^R/a and Ca = p^s/a, 
pi, /ii, a being the density, viscosity and surface tension of the dispersed fluid. U s and R 
are the jet surface velocity and radius, respectively. Alternatively, the Reynolds number 
Re = We/Ca = p\U s d/ ' pi and Ca can be used to characterize the jet dynamics; this choice 
is particularly useful in microfluidics, where Re is usually moderate or small (laminar flows). 
Surface tension is the main agent sustaining wave propagation of disturbances along the jet 
(downstream coordinate z), and thus Ca is the key parameter controlling the jet dynamics 
in microfluidics when Re is small. 



The link connecting connective/ abso 



well documented by experiments [7|, \8 



ute instability to jetting/dripping, respectively, is 



10] . Thus, the rate at which the perturbations grow 



is a function of Re and Ca. Since the jet disperses fragments of fluid 1 into an immiscible 
jet 2 of density p2 and viscosity /i2, two additional fundamental parameters are the fluid 
density and viscosity ratios p = P2/P1 an d p = pil P\- When jetting is produced by flow 
focusing 171]. one has a stream of fluid 2 forced through an orifice of diameter D which 
focuses the jet of fluid 1 (see Figure 2, inset). In flow focusing, we assume that R <C D, the 
most usual case. The relationship between the axial velocity U z of fluid 2 near the axis of 
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FIG. 1: (Color online) Different coflowing capillary jets into a much thicker silicone oil (which is in 
turn focused by an outer gas flow), from Ref (a) an air jet (/i = 2.06 x 10 4 ); (b) a mercury jet 
(fi = 2.42 x 10 2 ); (c) an ink solution jet (/i = 22.2). The lines guide the eye along the jet interface. 
Dash-dotted lines approximately delineate the positions of the centers of the cross section area, 
not in a straight line owing to a slight bend of the outer flow. Observe that the overall jet slope 
increases as the viscosity ratio /i decreases. 

the exit orifice, U z = U2, and the flow rate of fluid 2 through the orifice, Q2, depends on the 
Reynolds number Re^ = P2Q2/ {.P2D) (see Figure [2]). For Rei> — > 00, one asymptotically 



has|l8p U2 = O.5U2, where U2 = 4:Q2/(tiD 2 ). However, for Re^ 1, one approximately has 
U2 — 2t7 2 - For moderate to high Re/), a useful approximation is U2 ~ U 2 (see Figure [2] for 
Ken = 2200). Besides, when viscous effects dominate in the flow of fluid 1, the jet dynamics 
is thus determined by the parameters {Ca, p, fi} only, since Re disappears from the analysis. 
In the limit Re — > 00, Re is out from the analysis as well, and viscous effects are confined to 
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FIG. 2: Axial velocity profile u z = U z /U2 as a function of the radial coordinate r (made non- 
dimensional with D/2) for an incompressible fluid at a round orifice in an infinite plane thin wall 
as a function of the radial coordinate for various values. The flow is forced through the orifice 
by a pressure difference across the thin wall. In particular, the case Re^j = 2200 has been obtained 
using numerical simulation (Volumes of Fluid). The inset shows the basic flow configuration. 

boundary layers at the interface [3] . Since in this case the jet radius R is implicitly given 

by the equation p2^| = 2er/i? + p\U\, the jet dynamics becomes governed by {We, p, p, U}, 

where U = U 2 /Ui and U\ = Qi/(ttR 2 ) (Q 1 is the issued liquid flow rate of fluid 1). These 

two limiting cases allow for analytical treatment, and will be discussed in this work. 

Naturally, a jet is characterized by its slenderness and the applicability of the usual 

slender flow approximations. To study the linear stability of the system, we interrogate 

the jet response to small perturbation amplitudes £ <C R with a given axial wavelength 

A (see figure [3]). As long as the jet curvature is very approximately i? _1 (2), and the jet 

dR 

radius variations along axial distances A, of the order of O(RX) (where R = -3—), are small 

(Jj/O 

compared to R, one can always choose amplitudes such that O(RX) C 0(() < O(R), which 
justifies the classical cylindrical approximation. Under these assumptions, we investigate 
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FIG. 3: Sketch of a perturbed jet portion. 



the linear dynamics of cylindrical steady capillary jets in two limits: (i) creeping flow limit 
(Re — > 0), and (ii) inertia dominated limit (large Re). In the first limit, our exploration will 
focus on some practical scenarios of jetting at very small scales very recently reported in the 



literature 



10 



12 



20 



21 



22 



23]. 



Consider a thin cylindrical jet of fluid 1, moving in a co-flowing fluid 2. The general 
capillary flow category claimed in this work, for which the flow rate Qi could be made 
arbitrarily small (i.e., jetting would always take place independently of the issued flow rate), 
requires either making the jet radius R very small or allowing for recirculation (part of the 
jet inner flow should flow upstream) for continuity. Asymptotically, one may as well have a 
static jet core (U\ = 0) surrounded by a co-flowing fluid with velocity U2, which demands a 
negligible viscous diffusion from the interface. In this work, we will study these configurations 
which allow analytic dispersion relations. Thus, the occurrence of UJ is not restricted to 
either small or large Re values: in fact, the only necessary condition is that the jet surface 
velocity should be large enough to convect all perturbations downstream. The limit cases 
of large and small Re will be discussed in some detail in this work because they allow for 
analytical treatment and for their importance and illustrative power. In this sense, given the 
utter importance of controlled micron- and nano-sized droplet generation, we aim to provide 
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a global understanding of mechanisms supporting UJ, to guide future fluid disperser designs 
of special relevance in chemical engineering, combustion and energy efficiency, transport, 
food processing, spraying, biochemistry, pharmacy, biomedicine, environmental engineering, 
among others. 

A steady capillary jet of a fluid surrounded by an immiscible continuum fluid phase (fig- 
ured]), an intrinsically unstable state [yj], is locally stable (jetting is possible) whenever the 
so called marginal stability velocity v*_ relative to an observer is positive, so that all per- 
turbations are convected downstream. From now on, whenever we assume a local analysis 
applicable, velocities made non-dimensional with the local surface velocity U s . Since drip- 
ping is axisymmetric, here we consider axisymmetric perturbations only. Assuming small 
perturbations as superposition of waves proportional to exp[i(kz —ut)], where wave number 
k = k r + iki and frequency uj = uj t + iuji (made non-dimensional with R~ l and U s /R, respec- 
tively) are complex numbers, the jet linear dynamics is governed by the dispersion relation 
between k and u. Central to our analysis is the fact that steady capillary jets are unstable 
states [13I. \i\ . that is, they exhibit wave number ranges k such that perturbations grow in 
time as exp^i] leading to break up. A growing disturbance usually spreads along the jet 
bounded by two fronts moving with velocities v I and v* , the extremal values of the envelope 

nnnn 

velocities v = uoi/ki [13|, [14|, [15|, |16j. The extremal values or marginal stability velocities v* 



should satisfy 

v* = uj*/k* = duJi/dki\ k=k , , duji/dk r \ k=k * = 0. (1) 

Thus, the fate of the jet at a source-fixed station is determined by the minimum marginal 
stability velocity v*_. If v*_ > 0, small perturbations are convected downstream for all wave 
numbers (convective instability, or local stability), while if v*_ < 0, some wave number ranges 
will grow locally without bound (absolute instability). The aim of this work is to report 
a special class of parametric conditions of capillary jetting for which the marginal stability 
velocity v*_ (minimum front propagation velocity) keeps always positive for vanishing dis- 
persed flow rates of fluid 1. We designate this rather singular flow condition "unconditional 
jetting" (UJ). This means that the capillary jet is convectively unstable, or locally stable, 
and does not undergo a jetting-dripping transition as the issued flow rate Qi vanishes. The 
technological relevance of this class of flows can be understood as follows: picture a steady 
capillary jet flowing down from a slightly opened tap. If u* were always positive, one could 
slowly turn off the tap without transition to dripping. Eventually, when the tap is turned 



off, the jet would thin down to the continuum limit without transition to dripping. For 
extremely small flow rates, nearly monodisperse droplets of any imaginably small size would 
be produced upon Rayleigh breakup at a highly controllable rate. Although this behavior is 
indeed rather unusual for laminar jets from taps, it is however a real occurrence in co-flowing 
jets for a certain rather ample ratios of continuous-dispersed fluid densities and viscosities |4[] 
p and p, or when the continuous phase 2 co-flows with the jet at a velocity larger than a 

. nn 

critical velocity [12J, [19( u* 2 . We will see that the fundamental physical requirement for UJ is 
to have a jet surface velocity above a critical one. 

II. THE CASE OF NEGLIGIBLE INERTIA: ULTRA-THIN JETTING 



First, we will consider the case where velocities of both fluids are equal to the surface 
velocity (flat velocity profiles). Here, the jet's linear dynamics is governed by the following 



dispersion relation 



24j: 



iCaiu — k) 



N(k, uj, Re, p, p) 



2(1-/1) 



+ (k 2 - 1; 



D{k, uj, Re, p, p) 

For convenience, we define "viscous" wave numbers for both fluids as: 



0. 



(2) 



k\ = k 2 — i Re ■ (u — k) , k\ = k 2 — ip p 1 Re • (uj — k), 

Using these definitions, functions N and D are expressed as: 

N = 2kphk 2 [K Q {k 2 )h{k 1 )k 1 + I {k 1 )K 1 {k 2 )k 2 \ 
+k [k 2 (p - 1) - kj + pk 2 2 ] 2 h{k)h{k x )K Q {k)K x {k 2 ) 
+Ak 3 k 1 k 2 (^ - lfhik^hi^Koik^K^k) 
-hhik^Koih) { [k A + k\k\ + k 2 (k 2 - k 2 2 )] ph(k)K (k) 
+ [k{ + k\l - 2p) 2 - 2k 2 k 2 (p - 1)] I (k)K 1 (k)} 
-kJoi^K^) { [k\p - 2) 2 + 2k 2 k\p{p - 1) + p 2 k\] h{k)K (k) 

+ [k 2 (k 2 - hi) + k\{k 2 + k\)\ pI Q {k)K x {k)} 



(3) 



(4) 



D = k{[k 2 K Q {k 2 )K x {k) - kKoityK^h)] {kl - k 2 )h{k)I 1 {k 1 )+ 
ix \k x h{k x )h{k) - kh{k)h{k x )\ {k 2 2 - k 2 )K x (k)K x (k 2 )} 



(5) 



7 



Interestingly, a series expansion of N and D around Re = yields 

N = 0^1! at 2 ( k ,p)Re 2 + O (Re 3 ), D = ^ ~ ^ D 2 (k, /i)Re 2 + 0(Re 3 ), (6) 

where and D 2 (omitted, lengthy expressions) are independent of uj and p, as can be 
checked using Mathematica®. Moreover, in this limit the fluid velocities do not need to be 
uniform. Defining for convenience an average capillary number as 

C a = p 1/2 Ca = (pip 2 ) 1/2 U s /a, (7) 

equation [2] yields an explicit analytical expression for uj = u(k,Ca, p) providing a closed 
form for the frequency uj: 



uj = k + i(k 2 - l)Ca 



^N^p) + 2(/i _i/ 2 _ ^ 



(8) 



D 2 (k,p) 

This equation, together with conditions [fl provides v*, to* and k* for a given set {Ca, p}. 
Since Ca does not depend on the jet diameter d, a particularly strong practical implication 
follows: if one can find a parametrical region {Ca, p} where v* is always positive, it is so for 
any value of the jet diameter, no matter how small it can be. 

Figure H] shows a plot of the loci Ca = Ca*(p) where v* = 0. Fifteen orders of magnitude 
in p are explored, showing a small dependency of the critical Ca on p, which supports our 
definition choice for a relevant capillary number in the creeping flow limit (jet radii R <^ 
Pi/(piU s )), which incorporates both inner and outer fluid viscosities. This curve splits the 
{p, Ca} plane in two halves: above(below) this curve, v*_ is always positive(negative) and the 
jet flow is always supercritical(subcritical) independently of the jet diameter (supercritical: 
convective velocity always overcomes the upstream spreading of perturbations). Thus, if the 
velocity profile of both fluids is homogeneous in and out of the jet, supercritical jets of any 
imaginably small diameter could be produced for a co-flow speed u s larger than a critical 
velocity u* = Ca* (U* = oCa /(p 1 p 2 ) 1/2 )- 

Experimental results of other authors are compared with theory in Figure HI Anna 
et al.|20|] used a planar flow focusing device where they dispersed water (p± = 1 mPa-s, 
pi — 1 kg in silicone oil (p 2 = 5 mPa-s, p 2 ~ 0.9 kg L" 1 ). Their experiments show 



(Anna et al.|20J, Figures 3e,k,q) that jetting was found for values of the focusing oil flow 
rate Q 2 = 4.2 pLs~ l and above. Since their Reynolds number at the orifice was about 
6, a calculation of the oil velocity at the orifice axis yields about U 2 = 2Q Q /{hD) = 1.65 
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FIG. 4: Critical Capillary number as a function of the viscosity ratio fi. An asymptotic fit for // <C 1 
is provided: Ca* — > 0.1[ln(//)] 0,5465 . Also given: comparison with published experimental results 
were conditions for observed dripping/jetting transition are recorded, assuming a homogeneous 
velocity profile of both fluids in and out the jet. Unless otherwise stated, the estimated errors 
associated to data extraction from published plots are about ±20%. Data taken from (symbols, 



authors, data source figures): Q A nna et al. 



al. 



25|], Fig. 4; o Guillot et al. 



Fig. 3; □ Gahan-Calvo et al. 



20J , Fig 3; A Anna and Mayer 



231 ]. Fig. 5a; x Utada et 



ldl . Figs. 5(al,a2), and Herrada, Gahan-Calvo & Guillot 2008 



241 ]. Fig. 4. (liquid surface tension was 



app roximately 50 mN/m; 



liquid velocity at the orifice entrance can be calculated approximately 



26J] as U = kAP ■ D/fi 2 , 



where AP is the pressure drop through the orifice, and k is about 0.5). In particular, for Anna et 



al. 



201 ]. jjL\ = 0.001 Pa-s, /U2 = 0.006 Pa-s, 0.005< a <0.01 N/m (personal communication from the 



authors), 1.4< Q 2 <4.2 fjL/s or 0.27< U 2 <0.83 m/s, which gives 0.066< Ca <0.404 in this case. 



m/s, where h = 117 /im and D = 43.5 fim are the orifice depth and width, respectively (their 
orifice length was about L = 120 /im from their pictures, and thus a parabolic velocity profile 
should have developed). In accord with our predictions, they found jetting to occur above 
the indicated oil flow rate independently of the oil-water flow rate ratio, i.e. independently 
of how thin the jet was. Their corresponding threshold Ca = (ni^Y^^/o- is about 0.169 



(with errors associated to surface tension 



and indetermination between Q Q = 1.4/xL/s 
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and 4.2 /iL/s). Besides, using their same planar flow focusing device, Anna and Mayer [23] 
recently reported jetting independently of the focused flow rate beyond a capillary number 
Ca = 0.144 (worked out from their disclosed data), for an aqueous solution focused by oil 
with fi = 40, in the absence of surfactants. They reported transition from dripping to jetting 
for outer-to-inner flow rates ratios as large as 300, corresponding to jet diameters as small 
as about 3 /im. When surfactants are present, their results cannot be compared owing to 
non-linear dynamic surface tension effects beyond the critical micelle concentration c.m.c. 
(as they declare) at the jet. 

Moreover, Guillot et al. [l(| have reported an extensive and very valuable series of ex- 
periments of a liquid jet flowing coaxially in another immiscible liquid inside a cylindrical 
channel. When the jet to channel diameter ratio becomes very small, their measurements 
can be compared to our predictions. Their results show the dripping to jetting transition 
for a viscosity ratio// = 4.3 at Ca = 0.18 and 0.22 (calculated from their published data 
in their figures al and a2, respectively, for their higher available outer-to-inner flow rates 
ratios). We do not make use of their results with surfactants for the same reasons above 



given upon results from Anna and Mayer 23(j. Finally, for very large outer-to-inner viscosity 



fluid ratio, /i = 4.7 x 10 5 (syrup-air), we have found 24j transition from bubbling to jetting 
at Ca ~ 0.18. All these experimental findings, plotted in Figure IU provide full support to 
our prediction for homogeneous (flat) velocity profiles. 

Furthermore, we will consider a limit situation where the jet tappers from a finitely 
sized source (e.g. a capillary tube of inner radius R Q ), forming an arbitrarily thin spout of 
fluid 1 (figure EJ). For ease of understanding in this case, where we need to consider axial 
dependency of jet radius and velocities, we turn to dimensional expressions. Here, the local 
basic velocity profiles of fluids for negligible inertia, satisfying stress balance at the jet's 
surface, are given by: 

2" 



U 1 {r-z)=U s {z) + i 



4/t! 



1 



R{z) 



^H^)-^log(^) (9) 

To obtain these approximate equations, we have assumed that (i) the jet is slender (R 
is sufficiently small) and transversal velocities are neglected, and (ii) the outer pressure 
becomes negligible compared to the inner jet pressure p\ ~ a/R(z) as the spout radius 
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FIG. 5: Sketch of a jet tapering from a finite sized source in the form of an arbitrarily thin spout. 
In this configuration, the fluid velocities cannot be homogeneous for conservation of mass. 

becomes very small. Now, mass continuity for a vanishing issued flow rate of fluid 1, Q\ ~ 0, 
yields: 

v M = (10) 

Substituting this value of the surface velocity in the expression [7] of the average capillary 
number, one has that the jet would be locally stable if its local slope satisfies: 

• 8Ca* . . 

-R<- UT (11) 

This limiting slope is only a function of the viscosity ratio /x, plotted in figure [6j To obtain 
function R* = R*(/j,), we use the fundamental fact that the critical capillary number is 
independent of the velocity profile shape when Re <C 1, plugging in the values of the critical 
capillary number of figure HI This reveals that the jet slope is limited by the inverse of the 
viscosity ratio fi only, once the flow of both fluids is dominated by viscosity. It is worth 
noting that viscosity ratios smaller than about 0.01 would always provide, for Re < 1, a 
local stability of the tapering meniscus if the flow configuration allows co-flowing speeds U s 
larger than aCa* / (^i^) 1 ^ 2 , once jet is issued. 

The latter result, asymptotically valid for R small, is of fundamental importance to at 
least qualitatively explain experimental observations. In fact, if one observes figure [H the 
overall spout slope decreases as the viscosity ratio \i increases, according to equation (TTTT) . 
Moreover, it also explains why gas spouts are extremely difficult to achieve in coflowing 
liquids, since the local requirement ( TTTT) becomes very difficult to fulfill unless a very small 
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FIG. 6: Plot of the limiting slope of the jet as a function of the viscosity ratio \i. 

gas source is used. Nevertheless, upstream of the tapering thin spout, when the meniscus 
slope becomes of the order unity, the local Reynolds number may not be necessarily small and 
the limiting conditions for local stability deviate significantly from the above requirements. 
In fact, the critical capillary number decreases when Re increases (see Montanero & Ganan- 



Calvo29|, figure 5), and therefore the critical slope may increase, becoming of the order 
unity, in accord with real configurations like the ones shown in figure [D Notwithstanding 
this, again, requirement ffTTl) gives for the first time at least a qualitative explanation to 
observations on the slopes of tapering spouts and the difficulty to obtain stable gas spouts. 



A remark on the experiments by Courrech du Pont and Eggers|22j is here necessary: these 
authors report a "universal" critical value of the slope R = 0.47 ±0.06 of a tip singularity in 
a viscous combined withdrawal of air through a small orifice by means of a viscous liquid, a 
silicone oil with /i 2 = 30 or 60Pa- s. This remarkable result, from a large set of experiments, 
points to the existence of a locally self-similar conical region which eventually should tapper 
into an arbitrarily thin gas spout, where condition [TT1 should hold once the gas flow becomes 
viscosity-dominated in the spout. That self-similar region (we do not deal with this problem 
here) should be a function of the viscosity ratio fi only, although the structure of that 
locally self-similar flow is not known yet. However, our theory predicts critical slopes (about 
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6 x 10 4 ) much smaller than the one reported in 22], should a sustained (steady) thin 



spout be formed. Here we propose that the flow observed in [22j does not actually yield 
a steady extremely thin jet downstream of the observed region, within the suction orifice, 
for the experimental condition tested. In fact, these authors report that their setup "fails 
to produce a thin spout" , and that visible bubbles appear downstream of the orifice when 
entrainment takes place. The entrained spout would only be locally stable down to any 
imaginably small scale, leading to bubble sizes comparable to the spout diameter, if they 
had a local velocity at the orifice larger than 0.312 m/s or 0.441 m/s for their experimental 
conditions (fi2 = 60 or 30 Pa- s, respectively, with a = 0.0213 N/m and Ca* = 0.139). 
Since their velocities were below 0.0254 m/s (velocity at the axis corresponding to a flow 
rate of 0.01 ml/s through an orifice of 1 mm for viscosity dominated flow) they could not 
form a sustained extremely thin spout according to our prediction. However, the molecular 



mean free path for air enters into play below 1 micrometer 22j, and our theory becomes 
questionable below that scale. 



Finally, the elegant mathematical solution for a steady thin spout of Zhang 23] would be 
stable if her local convective velocities at the spout were everywhere larger than the critical 
ones here predicted. This condition would add an extra requirement in her analysis. 

The implication of present conclusions in microfluidics and, in particular, in the field of 
device design for emulsification are highly attractive for technological application. 

III. A CASE WITH DOMINANT INERTIA. 

While the above given analysis is restricted to negligible inertia, in the limit of dominant 
inertia one can also find an interesting, analytically tractable asymptotic case of UJ 4j, i. e. 
where the issued flow rate Q\ can be made arbitrarily small as well. In this alternative case, 
viscous effects are confined to thin boundary layers at both sides of the jet surface, where 
bulk velocities U\ and Ui are nearly flat but different. If both layers develop simultaneously 
from the same station near the dispersed fluid source, while they are thin compared to the 



jet 
in 



diameter the surface velocity U s can be explicitly expressed as (see the arguments given 

TT ITT 1 + (P^ 1/3U (,o\ 
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where, again, U = U2/U1. In this situation, if boundary layers are sufficiently small compared 
to the jet radius, the dispersion relation for an infinite cylindrical capillary jet moving with 
bulk velocity U\ in a co-flowing fluid with bulk velocity U 2 is 

/ m/ , ^o(*0 / , \ / , rr ^oW pU 2 k(k 2 -l) 

( w -k)(u- ku s )j^ + p{u - ku s )(u - kU)j^r = ±- L , (13) 

where We 2 = p 2 U 2 d/(2a). Note that the Reynolds number is absent in this expression, 
consistently with the assumption of dominant inertia. Viscous effects are however important 
at the jet surface, making p fundamentally important through the jet surface velocity u s = 
u s (U,p,p) 7^ 1. If we seek conditions yielding small dispersed flow rates, i.e. U — > 00, 
equation (fl~3]) plus conditions [7J with vZ. — yields the critical value We 2 using expansions 
u = Uq + U~ 1 lui and k = ko + U^ki (see notejsoj). Some flow configurations, such as 
flow focusing, give rise to further constraints; for very low viscosities, equilibrium requires 
We 2 — > 2 for U ^00 (from equation p 2 U 2 = 2a / R+ piU± ). In this case, for a given density 
ratio p, supercritical conditions are found for continuous phase viscosity larger than a critical 
ratio over the dispersed phase viscosity, or when p is smaller than a critical ratio for a given 



p value, as shown in figure [7] (see also {4], Figure 3) where thirteen possible flow focusing 
combinations are plotted. Note that the points correspond to given fluid properties, not 
reflecting transitions: the fluid combinations in the region "unconditional jetting" would 
exhibit this behavior, while those outside that region would only show dripping/jetting 
transition for a non-zero flow rate Q±. In other words, for the points in the region of UJ 
in figure [7J surface velocities are always large enough to have v*_ > and perturbations are 
flushed downstream, even though the bulk fluid can be literally static. 

IV. CONCLUDING REMARKS 

Interestingly enough, upon consideration of Ganan-Calvo 4J equation (8), the limit U — > 
00 here considered is strictly valid independently of the jet electrification as well, since the 
electrical term (finite) is overcome by other terms proportional to U ^> 1 and U 2 . Obviously, 
viscous diffusion of momentum from their surface soon affect the entire cross section of these 
jets. From well known boundary layer analysis, the downstream axial length that must 
be traveled to achieve Si ~ d should satisfy (piL^p^U" 1 ) 1 ^ 2 ~ Si ~ d. Thus, using equation 
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FIG. 7: Parametrical region in the {/x, p} space where unconditional jetting is found for flow 
focusing. Thirteen fluid combinations of interest are shown. Interestingly, water focused by air can 
theoretically exhibit unconditional jetting for ambient temperatures above 320K. 



121) . the maximum "unconditionally stable" length is limited to 



i " (w)' /3 „,> 



The equation for the surface velocity ( 1121) is then only valid within a length L M from the 
jet source, but not downstream. After this length, one cannot use the dispersion relation 
[13] since the bulk velocities become inhomogeneous; in fact, the liquid bulk should become 
recirculating in order to accomplish the vanishing Q\ condition (see figure [8]). Naturally, 
other cases of UJ may exist when inertia and viscous forces are comparable [e.g. that 
sketched in Fig. 8 (B)], but their required numerical treatment puts them outside the scope 
of the present analysis. 

As a final remark, it is worth stressing the fundamental importance of the surface velocity 
U s in the present analysis, for both creeping and inertia-dominated flows: note that for the 
cases here studied, where analytical solutions can be found, the wave behavior with respect 
to an observer is strictly determined by the linearized kinematic condition at the interface: 

where V is the small transversal component of the surface velocity. This is the only equation 
where the surface velocity appears: interestingly, since the first order velocity profile does 
not enter in the analytical dispersion relations here studied jsi|, when the observer moves 
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FIG. 8: Illustrative sketch of the UJ case for dominant inertia. (A) Flow within length from 
the jet source, where the asymptotic analysis is valid but viscous momentum diffusion from the 
surface thickens the boundary layers; (B) Development of a recirculation pattern downstream of 
the length owing to viscous diffusion. 

with the surface, its velocity disappears from the analysis, and one recovers the classical 
dispersion relations for cylindrical capillary columns existing in the literature. 
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